Double vector quarkonia production in exclusive Higgs boson decays 



Kartvelishvili, V./'B Luchinsky A.V., 2 'Q and Novoselov A.A. 2 '0 

1 Lancaster University, Lancaster, UK 
^Institute for High Energy Physics, Protvino, Russia 

Partial decay widths and branching fractions are calculated for the exclusive decays of the Stan- 
dard Model Higgs boson into a pair of vector quarkonium states, H — > J/tpJ/tfj, H — > TT,H — > 
J/ip(j>, H — > J/ipT, with relativistic corrections due to quark motion in mesons taken into account. 

PACS numbers: 12.38.Lg, 14.40.Gx, 14.80.Bn 



I. INTRODUCTION 

The necessary ingredient of the Standard Model (SM) 
is the Higgs boson — a scalar particle whose interactions 
are believed to generate masses of intermediate vector 
bosons and fundamental fermions. However, numerous 
attempts to discover this particle experimentally have so 
far only resulted in the limits on its possible mass Mh, 
114.4 GeV < M H < 182 GeV Q. 

Decay modes that can be used to detect this particle 
depend on Mjj- The coupling constant of the Higgs bo- 
son to another particle is known to be proportional to 
the mass of the latter, hence it is advantageous to search 
for Higgs boson decays into the heaviest of kinematically 
allowed particles. For instance, if Mh is above WW, ZZ 
production thresholds, H — > WW, ZZ decays (so called 
"gold-plated" modes) can be used. Observation of such 
reactions is one of the main goals of the LHC physics 
programme. 

However, the possibility of Mh < 2Mw cannot be ex- 
cluded. In this case bb and t + t~ modes are dominant, 
but large backgrounds makes these decays harder to ob- 
serve. In search for an exclusive Higgs decay mode with 
a good signature, in 0, [H, H[ decays of the Higgs boson 
into a pair of heavy quarkonia (H — > J/ipJ/ip, TT, J/tpT 
etc.) were considered. The widths of such decays were 
found to be small, but their good signatures and high at- 
tainable mass resolutions could make these decays useful 
in certain circumstances. 

The calculations in 0, H, H| were performed in the so- 
called (^-approximation, i.e. they did not take into ac- 
count relativistic corrections caused by the internal mo- 
tion of quarks in the vector mesons. Since then, how- 
ever, several theoretical and experimental studies of dou- 
ble charmonium production in exclusive processes have 
shown that the relativistic corrections can significantly 
alter the rates of these processes (see [!, [ll 0, H, H| for 
more details). In particular, the cross section of the re- 
action e + e~ — > JhpVc is increased by about an order 
of magnitude [1, E| 0|- Another example is the decay 
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Xb — * J/ipJ/ip'- in 0] it was shown that by taking the in- 
ternal quark motion into account, the width of this decay 
increases by a factor of 3. 

In this paper we study the influence of relativistic cor- 
rections on exclusive decays of the SM Higgs boson into 
a pair of vector quarkonia. The rest of the paper is orga- 
nized as follows: in the next section we briefly describe 
the formalism used, and present the distribution ampli- 
tudes for various vector quarkonium states. In section [TTT1 
analytical expressions are given for H — * ViV 2 decay am- 
plitudes corresponding to leading contributing subpro- 
cesses. In section ITVl we present numerical results for the 
Higgs decays into ViV 2 = J/ipJ/4>, T(15)T(15), J/># 
and J/ipY(lS) final states, and compare them to the 
results obtained within the ^-approximation. Our con- 
clusions are given in the final section. 

II. DISTRIBUTION AMPLITUDES 

Consider double quarkonia production in exclusive 
Higgs boson decays 

H -» V 1 { Pl ,X 1 )V 2 (p 2 ,X 2 ) 

where pi, 2 and A1.2 are momenta and helicities of the two 
vector mesons, respectively. 

In what follows, we restrict ourselves to the leading 
twist contribution. It can be shown that the final vector 
mesons in these decays are mainly longitudinally polar- 
ized, i.e Ai_2 = 0. Indeed, consider the explicit form of 
the polarization vector of a vector meson with energy E 
and momentum p. For a longitudinally polarized vector 
meson one has 

where My and Mh are masses of the vector meson and 
the Higgs boson, respectively. On the other hand, for a 
transversely polarized vector meson, 

e( A = ± i) = {o,-L, ± -L,o}~i. 

This vector is suppressed by a small factor ~ 
O (My /Mh) in comparison with the expression (p}, so in 
the following we will only consider the case \\ = \ 2 = 0. 
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The transition of a quark-antiquark pair into a lon- 
gitudinally polarized vector meson is described by the 
expression 



V(p,\ = 0) <f a (z)4(-z) 



f 5 ij 
1 



where a (/3) and i (j) are spinor and color indices of 
the quark (antiquark), respectively, while x is the quark 
momentum fraction with respect to the meson momen- 
tum. The constant / can be determined from the lcptonic 
width of the vector meson: 
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where a is the fine structure constant, e q is the quark 
charge (2/3 for J/ip, —1/3 for 4> and T), and M is the 
mass of quarkonium. 

A distribution amplitude ip(x), describing the internal 
motion of the heavy quark-antiquark pair inside their 
bound state, can generally be written as a series of 
Gegenbauer polynomials C^J 2 [Io| : 



6xx 



1+ Y, a n (»)C 3 J 2 (x-x) 
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where we have introduced a simplifying notation x = 
1 — x. The dependence of this distribution on energy 
scale pi is described by the QCD evolution of the moments 
a n (n): 



Mm) = ( J a n (yo), 



v a s (/io), 

with the anomalous dimensions 7„ defined by 
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where n/ is the number of active quark flavours. These 
anomalous dimensions are positive, so in the limit fi — > oo 
(or, equivalently, for light mesons) moments (0 tend to 
zero. As a result, in this limit distribution amplitude ^ 
tends to its asymptotic form tp a (x) — 6xx. We use this 
form of distribution amplitude for the <fi meson. 

However, for heavier quarkonia J/ip and T and for the 
scales /i = Mh ~ 100 -r 250 GeV, the asymptotic limit 
is not reached yet. Following [11(, for J/ip we use the 
distribution amplitude 



tpj/^(x,m c ) = c{(3) (l - (x - x) 2 ) x 

/ i 

^ 1 — [x — xy 



(8) 



with P 3.8 and the factor c(/3) fixed from the normal- 
ization condition 



ip(x)dx = 1 



o 



Note that this function is fairly close to a simpler expres- 
sion [H El 



(9) 



where ~ —3 is the intercept of the Regge trajectory 
corresponding to charmonium. For the wave function 
of T we use a parametrization similar to ^ with the 
intercept ay = — 9 The typical scale fi for these 

distribution amplitudes is the corresponding quark mass. 
Using formulae (g]), ([5]) and ©, their evolution to \i — 
Mh can be easily calculated. 



III. H 



Typical diagrams that give leading contributions to 
H — > V1V2 decays are shown in fig. [1] 

For the diagram shown in fig. [lp,, only Z will contribute 
as the intermediate virtual bosons, as the final mesons 
should be neutral and colourless. On the other hand, in 
this subprocess the two final mesons can be composed 
from different type of quarks. The corresponding ampli- 
tude is 



M {ZZr) = 



sin2?9 



w 



D Z {M 2 ) D Z (M. 



^M 2 H MiQ) 



Here $ w is the Weinberg angle, D Z (M) = M 2 - M% + 
iMT ' z is the inverse propagator of the virtual Z-boson, 
and constants ai,2 are defined as 



±- 



2e g sin d w 



with + (— ) sign used for c (s, b) quarks. One can see 
that this amplitude is expressed through experimentally 
observable quantities only, and hence does not depend 
on distribution functions (p(x). Note also that this am- 
plitude increases with increasing Higgs boson mass. 

The situation is rather different for the subprocess de- 
scribed by the diagram in fig. [1}d. In this case, either 
Z, W or gluons can be used as virtual particles. Their 
momenta can be calculated as 



xpi + yp2, 



q 2 = xpx + yp 2 , 



where p\^. are the momenta of final vector mesons, x and 
y are momentum fractions of quarks inside V\ and V2, 
respectively, while x = 1 — x and y = 1 — y . So, neglecting 
the masses of the quarkonia compared to Myy,z and Mh , 
the virtualities of the intermediate bosons are 



2xy{piP2) = xyM Hl 



<1~2 



1 ~ xyM 2 H . 
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FIG. 1: Typical diagrams describing H — * V1V2 decays. 



Clearly, in contrast with the previous diagram, these am- 
plitudes depend on the distribution functions (pi l2 (x). In 
addition, they do not increase with the increasing Higgs 
boson mass. 

In the case of ZZ decay, the quarks in mesons should 
be the same. The amplitude is equal to 



M {zz) 



1 3 



sin 2"dw 



b q M 2 H hf 2 M z x 



dx dy 



D z {xyM 2 H ) D z {xyM H 



(11) 



where 



=F 2e q sin 2 d w + Ae 2 sin 4 iV 



and the upper (lower) sign corresponds to mesons built 
from up- (down-) type quarks. In the limit of large Mh 
this amplitude remains constant. In the region Mh ~ 
2Mz, however, there is a noticeable peak, because in this 
case, due to the specific kinematics of the quarkonium 
formation, both intermediate Z-bosons are bound to be 
close to the mass shell. 

If the decay in fig. [TJj is mediated by a pair of W 
bosons, quarks in different mesons should be different, 
with one of the mesons built from up- type quarks, while 
the other — from down- type quarks. Here we consider 
J/ip(f> and YJ/ip final states. The amplitude of such a 
decay is equal to 



M (ww) 



1 e 3 cot dw , ,2 o 

-57 • 2 , v 12 My lhMz x 

24 sin v 



dx dy 



D w (xyM 2 H ) D w (xyM 2 H ) 



(12) 



where V12 is the respective CKM matrix element, while 
D W (M) = M 2 — Myy + iMT w is the inverse propagator 
of a virtual M^-boson. This amplitude also has a peak, 
this time at Mh ~ 2Mw- 

In the case of gg intermediate state, quark content of 
the vector mesons should again be the same. In the Stan- 
dard Model there is no H — > gg coupling at tree level, 
but it appears at higher orders when loop diagrams are 
taken into account. The main contribution to this effec- 
tive vertex comes for the top quark loop, for which we 
use the expression (see, e.g., [14l |) 

M(H -»• gg) = a °^ H \ I GfrGfr, 
D7T sm20wMz 

where the tensor G^ v is defined according to 



Strictly speaking, the value of the strong coupling con- 
stant in this vertex depends on gluon virtualities, but 
this dependence is only logarithmic and we do not take 
it into account. Using this effective vertex, we obtain the 
following expression for H — > gg — > V\ V2 decay ampli- 
tude: 



TV/to) = v^ etaniVa^A/g) 
27 M z 



2 x 



x / dxdy(p 1 (x)(p 2 (y) 
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xy yx 



(13) 



Clearly, this amplitude only depends on the Higgs boson 
mass through a s . 

There are other subprocesses contributing to the same 
double quarkonium final states, notably those contain- 
ing the tree-level vertex H — > qq shown in fig. [TJ;, first 
calculated in 0]. In such subprocesses, the two quarko- 
nia must again be the same, and the amplitude can be 
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FIG. 2: Amplitudes of different subprocesses contributing 
to the decay H — > J/ipJ/tp versus the Higgs boson mass. 
Solid (dashed) lines stand for internal quark motion taken 
into account (neglected). See text for further explanation. 



written as 



na s Ml hh 
9 M z M 



H 



dxdy 



xy 



Compared to other amplitudes considered above, in the 
Higgs mass range of interest this amplitude is suppressed 
by a small factor (My /Mh) 2 , so we will neglect it in the 
following. 



line). The amplitude of the subprocess shown in fig. QJ, 
(solid line labelled ZZr) remains unchanged when the 
internal quark motion is taken into account. As for the 
gluon-gluon subprocess (solid and dashed lines labelled 
gg), internal quark motion increases the corresponding 
amplitude roughly by a factor of 2. 

The main features of the above analysis also hold for 
other decays considered here. In fig [3] we present, from 
top to bottom, the Mh -dependence of the partial widths 
of the decays H -> J/ipJ/ip, H — > TT, H — > J/tpcj) 
and H — ► J/ipT, with the internal quark motion taken 
into account (solid lines) and neglected (dashed lines). 
The resonant peaks at Mh — 2Mz,Wi prominent in 
the ^-approximation in all these cases, are significantly 
smeared out once internal quark motion is accounted for. 
For Higgs boson masses away from these peaks, relativis- 
tic corrections increase the widths of the decays. 

Using the total width of the Higgs boson calculated fol- 
lowing [15l | , we have calculated the branching fractions for 
the four decays under consideration, from top to bottom, 
H -> J/ifjJ/ifj, H -> TT, H -> J/ip<t> and H -> J/tfjT. 
They are presented in figs. |4j It is clear that the branch- 
ing fractions vary strongly depending on the Higgs boson 
mass, and that the internal quark motion changes these 
fractions significantly. 



CONCLUSION 



IV. NUMERICAL RESULTS 

In order to obtain numerical results from eqs. (|10l — 
[T5)) . one needs the values for constants / and the dis- 
tribution amplitudes <f(x). The constants /, as deter- 
mined from the leptonic decay widths V — ► e + e~ using 
cq. ©, yield / = 200 MeV, 400 MeV and 700 MeV for 
(j), J/ip and T(15) mesons, respectively. Our choices for 
the distribution amplitudes of these three vector mesons 
were given by eqs. (jH [HI ©, respectively. After fixing 
these, calculations of partial decay widths for the decays 
H — > V1V2 are fairly straightforward. 

As an example, in fig. [2] we show M^-dependence 
of various amplitudes contributing to the decay H — > 
J/ipJ/ip. The amplitudes calculated with internal quark 
motion taken into account are shown with solid lines, 
while for dashed lines this motion is neglected. In 5- 
approximation, the amplitude of the subprocess H — > 
ZZ — > J/ifjJ/ijj has a prominent peak at Mh = 2Mz 
(dashed line labelled ZZ). The reason for this peak, 
as explained above, is that at this value of the Higgs 
mass both virtual Z-bosons in fig. [1J> are on mass shell. 
When the internal quark motion is taken into account, 
the quark momentum fractions x, y are allowed to de- 
viate from 1/2, and the Z-bosons are off mass shell. 
As a result, the peak is washed out, while away from 
the peak the amplitude becomes somewhat larger (solid 



We have calculated the partial decay widths and 
branching fractions for the exclusive decays of the Stan- 
dard Model Higgs boson into a pair of vector quarkonium 
states, H -> J/ipfaH -> J/ip, J/tj},H -> J/^T,H -> 
TT, with relativistic corrections due to quark motion 
in mesons taken into account, and compared our results 
with previous calculations performed in ^-approximation. 
We have found that the characteristic increases of the re- 
spective decay probabilities at Mh = 2Mz,w are still 
noticeable, but the relativistic corrections make them far 
less pronounced. For Mh away from that resonant area, 
relativistic corrections result in an increase of the branch- 
ing fractions by factors ranging from about 2 to 10 and 
above, depending on the specific decay mode and the 
value of the Higgs mass. However, the decay probabili- 
ties remain fairly small, and further studies are needed to 
identify possible areas where these decays may be useful 
experimentally. 
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FIG. 3: Partial decay widths for the exclusive decays H — > 
V1V2 versus the Higgs boson mass (in GeV) with internal 
quark motion taken into account (solid lines) and neglected 
(dashed lines): from top to bottom, the final states are 
J/ipJ/ip, TT, J/i>4> and J/V>T. 



FIG. 4: Branching fractions of the exclusive decays H — ► 
V1V2 versus the Higgs boson mass, with internal quark motion 
taken into account (solid lines) and neglected (dashed lines): 
from top to bottom, the final states are J/ipJ/ip, TT, J/ip4> 
and J/ipT. 



